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Abstract 

We generalise the classical Transition by Breaking of Analyticity 
for the class of Frenkel-Kontorova models studied by Aubry and oth- 
ers to non-zero Planck's constant and temperature. This analysis is 
based on the study of a renormalization operator for the case of irra- 
tional mean spacing using Feynman's functional integral approach. 
We show how existing classical results extend to the quantum regime. 
In particular we extend MacKay's renormalization approach for the 
classical statistical mechanics to deduce scaling of low frequency ef- 
fects and quantum effects. Our approach extends the phenomenon 
of hierarchical melting studied by Vallet, Schilling and Aubry to the 
quantum regime. 

Keywords: Transition by breaking of analyticity; renormaliza- 
tion; quantum scaling; specific heat. 



1 Introduction 

The Frenkel-Kontorova model (FK) is a one- dimensional lattice model 
exhibiting incommensurate structures. It is a system of elastically 
coupled particles in an external periodic potential (a discrete version 
of the sine-Gordon model) with Lagrangian 

L{X,±) ^'^l^ - v{Xn,Xn + l) \ , (1) 
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where 



V (x, x) = i (a;' - .t) ^ + P (a;' - a;) + cos {2-kx) , (2) 

with constants P ("pressure" or —P as tension) and u (amplitude 
of the onsite potential). We call successive pairs {x„,x„+i) bonds 
and V {xn, Xn+i) the (potential) energy of the bond. Notice that in 
Q and the particles' mass, the elastic coupling strength and the 
period of the onsite potential are all scaled to one. The FK model 
is a particular case of the broader class of (Generalised) Frenkel- 
Kontorova models (GFK) with Lagrangian still given by (0 but the 
bond action v is a generic function satisfying 

V (x + 1, a;' + l) = V (a;, x') 

In the space of parameters there are two important limits. In 
the integrable limit, the bond energy, v, depends only on (x' — x). 
For the FK model Q the integrable limit is attained at it = and 
its minimum energy configurations (i.e. a; G such that VM < 
N,Vm,n (x) = "^^^M V {xn, Xn+i) IS minimum for all variations of 
Xn with fixed xm and xn) are arrays of equally spaced particles with 
mean spacing 

Xn — Xm 
P ■= Imi — — 

-A/,JV-*oo N — M 

simply —P. In the anti-mtegrable limit 2] the onsite term dominates, 
which corresponds to it ^ 00 in (|5J. All particles are then in the 
minima of the potential and the mean spacing is the closest integer 
to —P, or any value in between the two if non-unique. 

An interesting set of codimension-2 critical points occur between 
these two regimes, often called Transition by Breaking of Analyticity 
(TBA): in the space of parameters (u, P), for each irrational p there 
is a curve P = Pp (it) of mean spacing p containing a critical value itc 
(there may be more than one Uc depending on the potential). The 
regime for it less than itc is called subcritical (or sliding phase) and 
above Uc supercritical (or pinned phase) For the FK model with 
fixed mean spacing p = 7"^, where 7 = (l + ^/E) /2 is the golden 
mean, the TBA is at the critical value Uc ~ 0.971 635 406 [B]. This 
is the case most often studied in the literature, since it is presumed 
to be the highest value of 11 at which a TBA occurs. 

For the case of mean spacing p — 7^^^, relabelling the bonds v 
appropriately (191 as r and v results in a Fibonacci sequence of rs 
and vs where each v-bond is always surrounded by r-bonds. The 
TBA point can be viewed as a fixed point of the renormalization 
operator that minimises the energy of the sum of successive v and 
r-bonds with suitably chosen space and energy scalings, a, J £ M. 
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respectively depending on the pair {v, r) (actually it is also necessary 
to subtract a constant and a quadratic coboundary but we will sup- 
press reference to these inessential terms) 1191 . The renormalization 
operator has a nontrivial fixed point^ (t*!^) with a ~ —1.414 836 
and J ~ 4.399 143 9. This fixed point corresponds to the critical 
Uc along the curve in {u, P) for p — 7~^, the transition point be- 
tween the subcritical and supercritical regimes. It has two unsta- 
ble directions: one along the curve of constant mean spacing, with 
eigenvalue S ~ 1.627 950 and the other transverse to this curve in 
the {u,P) plane (which we'll call the P direction), with eigenvalue 
V = -Jf-y ^ -2.681738 4 IT^I^ . 

Whereas the classical ground states of Frenkel-Kontorova mod- 
els have been extensively studied since the beginning of the 1980's 
121 El ISl 1221 (see also the review in chapter 1 of 1111 . and the 
book 1101 for several aspects of the FK model) , the extension to the 
quantum regime of the classical Transition by Breaking of Analytic- 
ity is still not fully understood. Most of the previous studies stem 
from the work of Borgonovi et al. |H1IH1, where the authors do a nu- 
merical study of FK model in the supercritical region u > Uc for the 
case of mean spacing 7~^. They introduce a 'quantum hull func- 
tion' for the expected positions the extension of Aubry's hull 
function (see 1111 . section 1.2) as well as the 'quantum g-function' 
(which reduces to sin (27r2;„) in the classical limit) and verify that by 
increasing h the quantum hull function becomes a smooth version of 
the classical hull function and that Qn tends to a sawtooth-like map. 
Similar results have been later obtained using various numerical or 
a combination of analytical and numerical methods 1^ 171 117l 1161 . 

In a recent numerical study Zhirov et al. 1261 claim to observe 
a 'quantum phase transition' in the FK model at a critical value of 
Planck's constant between 'sliding phonon gas' and a 'pinned instan- 
ton glass'. In fact we expect that due to KAM-type of arguments, 
for sufficiently irrational mean spacing the phonon energy band will 
survive for a small perturbation of the integrable limit u — 0, where 
the interactions between phonons of different wavenumber are small. 
At the anti-integrable limit, on the other hand, when the dominant 
interaction is the onsite periodic potential, the quantum spectrum 
consists of A*' degenerate sets of Bloch bands (A'^ being the total 
number of particles) whose width is due to tunneling or instanton 
effects between distinct minima of the onsite potential. As u is de- 
creased from the anti-integrable limit, the degeneracy between Bloch 
bands corresponding to distinct particles should be lifted, widening 
the Bloch bands until they merge into a unique phonon band for 
some non-zero value Uc (h) (possibly not a unique curve). 

In addition to the groundstates, it is physically significant to 

^This has now been proved by Koch 1181 by reformulation as a renormalization on 
continuous-time Hamiltonian systems and rigorous computer— assisted bounds. 
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study the effect of the TBA on the low temperature statistical me- 
chanics of FK models. This was done for the classical case in 1211 . 
The quantum statistical mechanics of FK models was considered in 
a series of papers by Giachetti and Tognetti, e.g. 1131 1141 . but we 
are not aware of any work on the effects of the TBA on the quantum 
statistical mechanics. 

The goal of this article is to study the transition above by extend- 
ing the minimum energy renormalization approach in 1191 to non-zero 
Planck's constant and temperature, which is done in section|5| This 
is done in a way similar to the classical non-zero temperature exten- 
sion performed in 1211 but requires first an extension of the classical 
case from ground states to time-periodic solutions. The strategy is 
to construct a renormalization operator, TZ, which reduces to the 
ground state operator, TZc^s, when h goes to zero, i.e. 7?.|^^q = TZcss- 
Then {u, P, . . . ,h = Q) is an invariant subspace for TZ which includes 
the critical fixed point of renormalization corresponding to the TBA 
by construction. 

In the main section |5| we start by defining a decimation proce- 
dure ® for pairs of "bond actions" by doing the trace over interme- 
diate particles corresponding to a partial trace over the kernel and 
the renormalization operator TZ by composing the decimation with 
the appropriate scalings. The quantum eigenvalue k is also intro- 
duced. Next, in subsection 12.11 we show that the classical ground 
state renormalization, 72.cgs3, is obtained as the limit of zero Planck's 
constant and frequency of quantum renormalization. In subsection 
12.21 the value of the quantum eigenvalue is determined by analysing 
the linearised or phonon problem. Finally in section r2.4l the results 
obtained are 'Wick rotated' to obtain the scaling laws for quantum 
thermal quantities, and an extension of the phenomenon of hierar- 
chical melting studied by Vallet, Schilling and Aubry |24l 1231 1^ to 
the quantum regime is proposed. 

2 Renormalization and scaling 

We begin by introducing the renormalization procedure for the action 
of time-periodic functions of prescribed period. Consider the class 
of models with the following formal sum for the action 

■,T). 

for time-periodic functions x„ of period T, where the bond action s 
is given by 

s{x,x';T) ^ 1^ _v(x,a;')|dt, (3) 
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and V is a GFK bond potential energy. The quantum renormalization 
will concern the trace of the kernel which is 



K' {T,h) = / P'xe^-^'"'^', 



where the notation D'x means that the integration is to be performed 
over the space of periodic paths x with period T and includes the 
integration over the endpoints dx (0) (i.e. x (T) = x (0) and D'x — 
da; (0) "Dx). Rescaling time t ^ t/T and defining — 2n/T the 
bond action can be rewritten as 

'^'"''t)^lfi (4) 

and now s acts on the space of period-one functions x, x' : R/Z — > R 
or loops. 

Now let X G R^ be a classical ground state and call a bond action 
s a r or an u as for the renormalization for classical ground states 
in )19| . For the case of mean spacing 7"^ the sequence of bond 
actions then forms an infinite Fibonacci sequence of r and v types 
of bonds where each v is always surrounded by rs. At this point 
one wants to eliminate all particles 2 from sequences of the form 
v{x,z;il) + t{z,x';Q) (for simplicity we use Q as an argument of 
action bonds instead of 2ti/Q, from now on). In order to do this 
define the following decimation operator © acting on pairs of bond 
actions (w,t) as 

{v®t){x,x'-Q) =-ifiln j D'^e^ (5) 

The decimated functional w © r still has the form of an action bond 
in the sense that it acts on pairs of loops {x, x'). The renormalization 
operator is defined as the composition of a decimation and scalings 
as 



7^ 



T {x, x'; n) 
V {x, x'; Q) 



J 

e 



{v © r) {x/a, x' /a; 0,/e) 
T {x/a, x' /a; fi/e) 



(6) 



which includes a scaling of frequencies, e, still undetermined at this 
point. The global scaling of bond actions m J je instead of just J 
because the renormalization now acts on actions instead of energies 
as in the classical ground state case. By performing the composition 
of and explicitly one can see that the natural scale factor 
K = J je for Planck's constant arises and TL can be seen as acting on 
the extended space of (r, w, K) for fi > as 

t(x,x'-^) = -iKfi/D'2e^^{"("''"'-"/"'+K-'-'/<«;n/^)} 
Tl:\ v{x,x';Q.) = ^T{x/a,x' /a;Q,/e) 
fi — nfi 

(7) 
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We interpret k as the eigenvalue of TZ in the direction of Planck's 
constant. 



2.1 Semiclassical approximation and ground 
state limit 

As a first step to connect the quantum renormalization @ with the 
ground state renormalization of |19| . one can look at the decimation 
operator @ in the stationary phase approximation for small fi: 

{v®t) {x,x;n) = -ifiln J V' z S [z - Zd. {x,x';n)) x 

xei{"(^'^;")+^(^'^''")}. (8) 

Here Zd. {x,x';il) is the classical path of period one satisfying the 
Euler-Lagrange equations 



Slv {x,z-Q,) +T{z,x';n)] 



5z 



= 0. (9) 



In the classical limit one is thus left with a dynamical problem Q 
and the decimated action JSJ can be rewritten as 

sta \v{x,z;fl) + T (z,x';Q,)] , (10) 

the sum of bond actions evaluated at z = z^i., given by which 
stationarises the sum v + t over the space of all loops {z : R/Z — > R}. 



2.1.1 Ground state limit 

The ground state decimation can now be taken as the limit f2 — > 
of or lion . To see this, notice that in this limit only the classical 
ground states contribute to the kernel and 

|,,,*(**,(.V„,,-»,...-.,„H..-.(.,.0), 

where (x, z,x') is a segment of a classical ground state and z(x,x') 
minimises the sum v + t. Taking the logarithm and multiplying by 
— ifi, this results in the decimation Q being 

(ver) (x,x';0) = lim -^min [v'"' (x,z) H-v'"' (z,x')l (11) 

(the factor 2tt/^1 is kept in the above expression simply to con- 
trol the divergence of decimation). Apart from the limit factor 
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limn^o —2n/Q, this is in fact the classical ground state decimation 
in I19| and so the renormalization © becomes 



7^ 



r (a;, x'\ 0) 
V [x, x \ 0) 



lim 

S2^0 



vW) (x/q.x'/q) 



-2f:rv(-) (x/a,x7a) 



27r^ 
lim --=-/<.cgs 



(x,x') 
(x, x') 



(12) 



where ©cgs and 7?.cgs are the decimation and renormalization opera- 
tors for the classical ground states, and e £ R is still undetermined 
at this point. The renormalization TL has therefore a fixed point (the 
TBA fixed point), with ground states of mean spacing p — 7~^, in 
the f2 = subspace consisting of 



r(a;,x-';0) 
V {x, x'; 0) 



lim — ^ 



{x,x') 
ix,x') 



where (v"^, v^) is the fixed point of the ground state renormalization 
operator 7?.cgs, because applying expression 11211 at the TBA fixed 
point results in 



7^ 



r (a::, x'; 0) 
V [x, x'; 0) 



= lim — — 



(x,x') 



T {x, x'\ 0) 

V {x,x'; 0) 



2.2 Renormalization for the phonon spectrum 

Close to the ground state, for small , the relevant contributions 
are approximately given by the normal modes or phonons. For a 
GFK model at irrational mean spacing ground state, p, the phonon 
spectrum includes zero in the subcritical regime, but the minimum 
frequency or phonon gap is positive above the critical point Uc 0- 
The phonon contribution to the trace of the kernel is given by doing 
a quadratic approximation which is (with Xn — x„ + ^„) 



where the kinetic term was generalised to a symmetric quadratic 
form in the velocities 



e'M ."" + 1) X (13) 



(as will become clear later on the renormalization operator intro- 
duces coupling between the velocities of neighbouring particles) , and 
u is also a symmetric quadratic form 
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The first term in 1131 is simply the classical ground state and 
corresponds to il'21 . so one wants to apply the decimation Q to the 
sum of the quadratic parts of bond actions of the form 

S {iJ (c,c) +-'^^ (^'^')] - H"' (^'^')] 

Because this sum is quadratic in the particle to eliminate, (, the 
corresponding functional integral can be calculated (for example 
by Gaussian integration of the Fourier transformed sum of action 
bonds |27| 'l and the semiclassical approximation is exact in this case. 
To first order in Q,^ the result is^ 
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Kfi{i&K"'(«^c)+-<-'(c.e')]~["'"'(«.o+..<-)(c,«')]} 



:e'^ 



i|^{ii^„-A)(j,e')_i,(.)(5,«')} 



where m'-^^ is a new symmetric quadratic form with components 
given by 

m l' = m[V + , " "7 , + ^ 



and the new form ft'-^-' has components 



(u) 

-(r) (u) 
Ml = U)^ 



(v) I (t) 
U22 ' Mil 



(t^) (r) 

~(t) M12 M12 /I p:\ 

"12 - (t;) (r) 

17(22 "I" 1 

(r)2 

-(^) _ (t) _ M12 
M22 — M22 



M22 ~r Mil 



With these new quadratic forms, the decimation becomes simply 
{v (Br) (^x,x';Q,) — (c.g.s. decimation) + (16) 
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There is possibly also a logarithmic term due to the integration measure which 
corresponds to a redefinition of the ground state energy. 
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where 'c.g.s decimation' is the decimation (IIH . For the renormal- 
ization (|SJ| one also needs the 'undecimated' bond actions which cor- 
respond to isolated bond actions of type r, so one should define also 



(17) 



For a GFK model Q at the start of renormalization, for either type r 
or V of bonds the quadratic forms are mi2 = 0, mn +77122 = 1, Ujj = 
v,jj (x, x') for J = 1, 2 and «i2 = — v,i2 (x, x'), where the subscripts in 
V denote differentiation with respect to the first and second variables 
and (x,x') is a segment of a classical ground state. Because u is then 
dependent on (x,x'), after iterating the above transformations 11411 . 
HSU and 1171 1 one ends up with a set of asymptotic quadratic forms 
depending on the ground state, m'^', fh!-^\, u'^', and v}-"\ which 



scale by factors'^ 



a e / J ~ 1.255 071 for the mass forms and 



J /a for the potential forms EHli i-e. 



m. 



J x,x' 
= (v) a?!p- - (t) 

x,x' x,x 

x,x' J/ x,x' 

x,x' J x,x' 



(18) 



Finally, using 114m8L and defining 

f[x,x\VL) ■- ^ [-v<^' (x,x') 
\ ( rfl 



+ 



v{x,x'-M) := g[-v(")(x,x') 



4^2 x,x' ; 



the renormalization (O at (r, ti, /i) becomes 



7^ 



r (a;, a;'; fi) 
V (x, X \ Q.) 

n 



f {x, x'; Q) 
V {x, x'; ^1) 



Thus, by including the scaling of frequencies e ~ 1.649 415 (see foot- 
note |2J in the renormalization ®, the point {f,v,0) becomes an 



^ In )19l the scaling is actually for the quantities an 



,("-!) 



„(") 



and d„ := m^j' with scaling constants oj ~ 1.255 071 for a„ and 
bn and /3/a for Cn and dn, so the first and third equations in 1181 are defined up to 
a quadratic coboundary. Here we use the scale factors for frequency e := ^ldJ jop- ~ 
1.649 415 and energy J = a/3 instead. The origin of is still a mystery. 
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approximate fixed point for small fl (actually a line of fixed points 
parametrised by Q). In particular the phonon spectrum is asymp- 
totically self-similar under scaling by e in the direction of Q,. 

Finally, note that fixing e also determines the eigenvalue in the 
fi-direction* as k = ^/e ~ 2.630 716. 



2.3 Scaling of the kernel 

The results above imply that the eff'ect of including both the fre- 
quencies and quantum directions, close to the TBA fixed point (for 
Au := u — Uc, AP := P — P^-i (y-c), and H small), is that the 
following asymptotic relation of bond actions (regarded as functions 
in parameter space) holds 

(r, v) {Au, AP, fi) ~ J (r, v) (SAu, t)AP, e^, kH) . (19) 

If K'j is the trace of the kernel for a chain of size Fm, the mth 
Fibonacci number, in the discretised form with J 'time steps' (such 
that K'f^^ ~ limj^oo K'jp^^), then 



K'j,p^ {Au, AP, n, h) ~ K'j,p^ ^ {5 Au, vAP,en,^fi)(^ 

Here the multiplying factor comes from the functional integration 
measure due to the change of coordinates (here with diagonal mass 

. (n) - (n-1) , _ (n)N 

components fi^ ' — -f mji ) 



j-i 

X, 



n(^:^^)(".^o - nnv 4.^i, 

n 71 j—0 

= nfjS) {'P'jin){en,^n), 

where x„ are the positions of the renormalized particles. 



2.4 Non-zero temperature scaling 

The quantum partition function at temperature can be easily ob- 
tained from the trace of the kernel by Wick rotation, i.e. doing 
27r/f2 = —ipti, where (5 — I/O. The quantum partition function is 
then 

Z{l3h,h) =K' {-il3h,h) =y D'a;e-'''^^'"''"''\ 

*This result was published in 1121 containing a mistake: k = Je instead of the correct 
value ft = J je. 
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where the Euclidean action for the FK model is 



Se = Se {x„,x„+i; ph) 



= E 



2 ipny 



and the integration is now over the rescaled Euchdean time ts G 
R/Z. Comparing with the Euclidean bond action Se can be 
written as 



Se [x,x';f3H) = 



— — s (x, x'\ Q.\ 



and has units of an energy instead of an action because of the multi- 
plying factor Q,/2ti. This gives another direction to which renormal- 
ization can be extended, the temperature direction G = Define 
then the operator ©b by 'Wick rotating' the operator © for the trace 
of the kernel @ , to obtain 

{vE®ErE)[x,x'-l3h) = -p-^\n J i)'^ e-''["«'"'"''^")+"^("'"'''"'')l , 



and the renormalization as 
TZe 



Te [x,x';Ph) 
Ve {x, x';/3h) 



Jte {x/a, x' /a; ef3h) 
Jve {x/a, x' /a; e/Sh) 



Therefore, close to the TBA fixed point the renormalization opera- 
tor leads to the corresponding relation to 11911 for Euclidean bond 
actions: 

{TE,VE){Au,AP,/3h,R)::^J{TE,VE) {SAu,r,AP,e-'f3h,Kh) , 

so under the renormalization TZe there is a fixed point at = 0, 
/3fi = oo, with an unstable eigenvalue of ree = ^7 ~ 4.339 143 9 in 
the temperature direction. Although this agrees with the scaling of 
for classical statistical mechanics, the result here includes the 
momentum contribution and corresponds to the quantum correction 
of the low temperature classical result, in the region O ^ in which 
the classical partition function is not valid (see 1111 . section 3.2). 
Denoting by Zf„^ (Ait, AP, 0, h) the quantum partition function for 
a chain of Fm particles at temperature Q, regarded as a function in 
parameter space, the renormalization picture leads to 



Zf^ {Au,AP,e,n)~ZF^_, {5Au,r)AP,JQ,Kn) 



The free energy / = — 01im„i-.oo In^F^/i^m therefore behaves like 
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T 



Figure 1: Example of the behaviour of Cp for the choice of k{a,b) — 
[c+fc+ (a, 6) + c'k' (a,6)]e"''(T^"i^) in equation (|^ with c = 1/2, c+ = 1, 
= 1/5 and A;± = 1 + cos27r Ij^ ± 

' V In J/ m K / 



(using 7 ^ = limm^oo -F,ti-i/-F,ti) 

/ (Au, AP, e, fi) ~ 4?/ ('5 A«, 77 AP, ^ e, fi) - — In (ff) , 

close to the TBA fixed point, and e — — d (//O) /96, the energy 
per particle, like 

e (Aii, AP, e, h) ~ — e (S Au, r) AP, JQ,K.h). 

This leads to the following scaling law for the specific heat per par- 
ticle at constant pressure Cp = de/dO: 

Cp ( Au, AP, e, fi) ~ 7~^cp ((5 Au, riAP,je,Kn). (20) 

In particular close to the TBA, with Au = AP — for small O and 
h, equation H2()^ suggests that in the (0,fi)-plane the specific heat 
shows a sequence of modulated ridges, which are invariant under 
scaling by [J' ,k): 

Cp (e,?i) ~ OTS^ (In e. In fi), (21) 

with fc (a + In J", fo + In k) = k{a,b), for all a and b (see figure Ql. 
These ridges thus constitute the extension to the {&, fi)-plane of the 
sequence of Schottky anomalies responsible for the phenomenon of 
hierarchical melting of the chain |24l 1231 125] , 
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